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INTRODUCTION 


This  report  addresses  some  theoretical  problems  associated  with  the 
railgun.  It  Is  comprised  of  five  sections.  In  the  first  section  a  fulLy 
relativistic  analysis  is  made  of  a  simple  prototype  of  the  railgun. 
Electromagnetic  field  criteria  are  obtained  for  the  realization  of  maximum 
propulsion  speed.  The  second  and  third  sections  address  an  equation  of  state  for 
a  strongly  coupled  pLasma.  Results  relevant  to  two-component  plasma  given  in  the 
second  section  are  generalized  to  a  multi-component  plasma  in  the  third 
section.  A  significant  result  of  this  work  is  exhibited  in  a  plot  of  PV/NkgT  vs. 
temperature  T  which  exhibits  variation  with  the  shell  structure  of  atoms  in  the 
plasma.  Application  is  made  specifically  to  an  aluminum  plasma. 

Studies  of  electrical  conductivity  in  a  strongly  coupled  plasma  are  reviewed 
in  the  fourth  section,  and  in  the  final  section,  an  analytic,  expression  is 
obtained  for  nonlinear  electrical  conductivity  in  strongly  coupled  plasmas. 
Figures  and  references  as  included  at  the  end  of  each  section. 


MS 


FUNDAMENTAL  RELATIVISTIC  SOLUTION  FOR  A  RAIL  CUN 


Summary 

J- 

A  fully  relativistic  analysis  is  made  of  the  dynamics  of  a  rail"gun  based 
on  three  assumptions:  (1)  Ohm's  lav  is  valid  in  the  rest  frame  of  the  plasma.* 
(2)  total  electron  momentum  is  transferred  to  the  projectile.*  and  (3)  motion 
of  the  projectile  is  constrained  to  one  direction.  With  these  assumptions,  a 
relativistic  equation  for  the  velocity  of  the  projectile  is  obtained  whose 
solution  monotonically  increases  to  one  of  two  values  depending  on  field 
strengths.  For  B  >  E,  the  maximum  velocity  is  cE/B  whereas  for  E  >  B  it  is  c 
where  c  is  the  speed  of  light,  and  E  and  B  are  applied  electric  and  magnetic 
fields,  respectively  (in  cgs).^ 


INTRODUCTION 


Attention  has  recently  been  given  to  dynamical  properties  of 
electromagnetic  propulsion.  In  the  rail-gun  device.^  a  non-conducting  pro¬ 
jectile  is  propelled  by  a  current-carrying  plasma  driven  by  the  Lorentz  force. 
Plasma  dynamics  is  more  difficult  due  to  ablation  of  the  projectile  and  for 
the  most  part  previous  studies  have  attempted  to  incorporate  this  effect. 

In  the  present  study,  we  return  to  a  more  elementary  configuration  for 
the  purpose  of  developing  a  fully  relativistic  study  of  this  problem.  Thus* 
for  example,  it  is  assumed  that  the  rest  mass  of  the  projectile  is  constant, 
and  that  total  electron  momentum  from  the  plasma  is  transferred  to  the  projec¬ 
tile.  Furthermore,  it  is  assumed  that  the  projectile  is  constrained  to  move 
in  one  direction.  Our  remaining  assumption  is  that  Ohm's  lav  is  valid  in  the 
rest  frame  of  the  projectile.^ 

With  these  assumptions  at  hand,  a  relativistic  equation  is  constructed 
for  the  projectile  velocity.  Solution  to  this  equation  reveals  tvo  asymptotic 
velocities  which  depends  on  initial  field  strengths.  Thus,  for  example,  for 
the  case  E  >  B.  the  velocity  is  c.  the  speed  of  light,  whereas  for  B  >  E,  the 
velocity  is  cE/B.  It  is  further  demonstrated  that  for  initial  velocities  less 
than  respective  asymptotic  values,  velocities  monotonically  approach  their 
respective  limiting  values.  For  the  case  B  >  E.  starting  velocities  greater 
than  cE/B  are  found  to  decay  to  this  asymptotic  value. 


ANALYSIS 


Our  starting  equation  is  Ohm's  law*  which  in  the  rest  frame  of  the 
projectile  (primed  coordinates)  is  written 

£’ (i) 
where  a  is  conductivity.  Transforming  back  to  the  lab  frame  (see  Figure  1)  we 
find 


Y(J  -  c  B  p)  =o-E 

X  X 

J  =o-  Y(E  -  B  B  )  (2) 

y  y  r  z 

J  -<r  Y(E  -  8  B  ) 
z  z  r  y 

2  -1/2 

where  j)fv/c  and  Y  =  (l-{)  )  .  The  charge  density  p.  for  a  charge-neutral 

plasma,  is  equal  to  2ero.  In  the  lab  frame  we  take 

1-4  (3) 

_  A 

B  =  Bz 

where  hatted  variables  denote  unit  vectors.  Inserting  these  values  into  (2) 
gives 

J=0’Y(E-jJB)y  (4) 

with  J  denoting  microscopic  electron  velocity  we  write 

J=qn<5>  (5) 
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Electron  charge  and  density  are  q  and  n  respectively.  We  further  recall  that 
the  density  transforms  as 


m 
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m 

Civ 


a  s  Yn' 

Combining  the  latter  four  equations  gives 


(6) 


<15  >  s  V  (E  -  p  B)  1 


(7) 


where 


U  = 


qn' 


(8) 


represents  mobility  in  the  rest  frame. 


Taking  the  average  of  the  Lorentz  force  on  electrons  we  obtain 


•4-< 


-  A. 


,n  >  =""*■<  Ti  >  B 
dt  py  c  *y 


(9) 


where  we  have  recalled  the  vector  property  given  by  (7).  We  assume  a  total 
transfer  of  electron  momentum  to  the  projectile  which  gives 


p  >  y  =  ““  P 

dt  *y  ~  dt  ~ 


(10) 


where  N  is  total  number  of  current-carrying  electrons  in  the  column.  The 
momentum  of  the  projectile,  P  •  is  given  by 


P  =MYv 


(11) 


where  M  is  the  mass  of  the  projectile. 


Combining  (7),  (9),  (10),  and  (11)  gives  the  desired  equation  of  motion 


dt 


YMv  =  -aL(  1  -X)  B 


(12) 


vhcre  a  is  Che  length  of  che  conducting  column.  I  is  the  current 


end 


al  =  NqyE 


_  xS. 

w  = 

“  B 


(13) 


(14) 


Note  in  particular  that  from  (13)  ve  may  write 


E  = 


A  o* 


(15) 


where  aA  represents  the  volume  of  the  conducting  column  in  the  rest  frame. 


Integrating  (12)  gives 


0  c  **v 


n  -  n 


where  n  is  the  dimensionless  time 


du 


[  1  -  (■U*)Z]3/2  (1  -u) 

u  c  J  w 


(15) 


n  f 


_  aiat 


Me' 


(16) 


and  v=v  at  n  .  From  (15)  we  see  that  T?*oo  at  the  singular  paints  u  =  c 
o  o 

and  u=w  which  represent  asymptotic  velocities.  It  will  be  shown  below  that 
these  asymptotic  velocities  are  approached  monotonically.  With  this  property 
we  may  conclude  that  for  zero  starting  velocities  maximum  values  are  given  by 


E  >  B 


v  =  c 

max 


(17) 


v  =  w 

max 


E  <  B 

A  sketch  of  these  findings  is  shown  in  Figure  2 


To  examine  the  monotonicity  of  v(t)  we  differentiate  (12)  to  obtain 


(18) 


We  conclude  chat  for  vie  and  v  i  w,  dv/dn^o.  Furthermore,  with  v=o  at  t=o. 
(18)  gives  the  starting  acceleration  (in  dimensional  form) 


dv  _  alB 

dt  o  ~  Me 


(19) 


Note  that  for  the  case  E  <  B,  an  initial  velocity  vq  >  v»  decays  to  v  =  w  as  is 
evident  from  (18).  Furthermore,  as  is  clear  from  (15).  asymptotic  values  (17) 
are  independent  of  initial  velocities. 

Characteristic  times  corresponding  to  the  maximum  velocities  (17)  are  as 
follows.  In  the  limit  w  >>  c.  (12)  gives  the  characteristic  time 


with  maximum  velocity 


hcL 

T1  =  alB 


v(1)=c 

max 


In  the  limit  w  <<  c.  (12)  has  the  solution  (with  v=o  at  t=o). 


v  =  w  £l  -  exp(-t/T2)J 


where 


and 


r  =-S-t 
2  B  1 


(2)  . 
v  =  w 
max 


We  note  that  although  'r2<<rl’  accelerati°Q3 


(20) 


(20a) 


(21) 


(22) 


(22a) 


(1)  (2) 

V  V 

max  max 

T1  =  r2 


are  the  same* 
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APPLICATION 


In  applying  Che  preceding  results  to  experimental  values,  first  we 
rewrite  w  in  practical  units.  Setting 


(19) 


where  V  is  the  potential  across  the  rails,  permits  (14)  to  be  rewritten 


cV 
w  =  ~r 
aB 

In  practical  units  this  expression  becomes 


(20) 


/J&SL  _ 

w  (  )  * 


■LQ.  V(KY) 


s  B(kgauss)a( cm) 


(21) 


5.7 


Typical  experimental  values  are  :  V  =  lkV.  a  -  1cm,  B  -  200  kgauss  which 
gives  w  -  5km/s.  This  value  agrees  in  order-of -magnitude  with  observed  max¬ 
imum  velocities. 

CONCLUSIONS 


We  have  examined  the  relativistic  solution  to  the  rail-gun  configuration. 
Incorporating  some  simplifying  assumptions  we  found  that  the  projectile  velo¬ 
city  goes  monotoni cally  to  the  minimum  of  the  two  velocities,  c  and  cE/3.  The 
asymptotic  value  c  corresponds  to  E  <  B  whereas  the  value  cE/B  corresponds  to 
the  limit  3  >  E.  It  should  be  emphasized  that  this  present  study  does  not 
take  into  account  thermodynamic  effects  such  as  momentum  imparted  to  the  pro¬ 
jectile  from  the  exploding  "fuse.” 
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ONE-COMPONENT  PLASMA  AND  LARGE  Z  LABORATORY  PLASMA 


Summary 

Stemming  from  an  expression  for  the  mean  two-particle 
potential  energy  of  a  two-component  laboratory  plasma  comprised 
of  ions  of  charge  Ze  and  electrons,  it  is  argued  that  for  suf¬ 
ficiently  large  Z,  thermodynamic  properties  of  such  plasmas  are 
the  same  as  those  relevant  to  Wigner's  one-component  plasma 
model.  Thus  the  following  equation  of  state  is  obtained  for  a 
laboratory  plasma  with  Z  >  5  and  r  >  1 . 


nkBT 


l  +  4  +  t  (ar  +  br1/4  +  cr-1/4) 


In  this  equation,  P  is  pressure,  T  is  temperature,  n  is  ion 

number  density  and  a,  b,  c,  d  are  known  numerical  constants. 

2  1 

The  plasma  parameter  F  =  (Ze)  /akgT  where  4ia  /3  represents  mean 
occupation  volume  per  ion. 


15 


Strongly-coupled  plasmas  play  a  role  in  recombination 
approaches  to  x-ray  lasing  [1],  inertial-confinement  fusion 
devices  [2],  the  interiors  of  certain  super-dense  stars  [3], 
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and  in  plasma-driven  rail-gun  devices  [4]. 


For  the  most  part  studies  addressing  strongly-coupled 


plasmas  [5,6,7]  have  employed  a  plasma  model  due  to  Wigner  [8] 


which  was  conceived  for  the  purposes  of  studying  phase  transi¬ 


tion  to  the  solid  state:  This  medium  is  called  "jellium",  or 


more  commonly,  a  "one-component  plasma",  which  often  carries 


the  abbreviation,  OCP .  An  OCP  is  comprised  of  ions  moving  in 


a  charge-neutralizing  uniform  negative  background. 


In  the  present  work  attention  is  directed  at  a  laboratory 


plasma  comprised  of  ions  of  charge  Ze  and  electrons  in  equilib¬ 


rium  at  a  given  temperature.  Examining  the  interaction  energy 


of  the  plasma  indicates  that  for  Z  >  5,  thermodynamic  properties 


such  as  internal  energy  and  equation  of  state  are  given  by  cor¬ 


responding  expressions  appropriate  to  an  OCP.  As  an  application 


of  this  finding,  the  equation  of  state  for  such  a  relatively 


high  Z  plasma  is  obtained  from  a  previously  constructed  expres¬ 


sion  for  the  Helmholtz  free  energy  [9] 


With  the  plasma  under  consideration  comprised  of  electrons 


and  ions  of  charge  eZ,  charge  neutrality  implies  the  constraint 


eV(ng  -  Z  n^)  =  0 


where  V  is  plasma  volume  and  n  and  n_  are  electron  and  ion 
r  e  L 


number  densities,  respectively.  The  aooroximate  eaualitv  in 


ms 


r.  r.  <.  * 


p**Kv.x  .x 


<•1 


(1)  derives  from  the  inherent  statistical  nature  of  a  plasma. 


For  a  two-component  plasma,  three  interactions  contribute  to 


the  mean  two-particle  potential  energy  and  we  write 


<V>  ,  (4»)1/3[nyv 


(nenZ>1/6Ze2  +  n3/3(Ze)2]  (2) 


With  the  constraint  (1)  we  may  set  ng  *  Z  nz  and  (2)  becomes 


<V>  =  e2(47rnz)1/3(Z1/3  -  Z7/6  +  Z2) 


Thus  in  the  limit 


Z2  »  Z1,17 


the  relation  (3)  reduces  to 


<V>  =  (4irn2)1/3(Ze)2 


which  we  recognize  to  be  the  interaction  potential  of  a  one- 


component  plasma  comprised  of  ions  of  charge  Ze. 

This  similarity  may  be  further  illustrated  through  the 
plasma  parameter  [10,11] 


2/3  <V>  <V> 


where  <E„>  is  mean  kinetic  energy  per  particle.  The  relevance 


of  y  Co  the  properties  of  a  plasma  is  evident  from  (6).  Namely, 


with  this  expression  we  may  conclude  that  a  plasma  is  strongly 


coupled  when  y  >  1  and  weakly  coupled  when  y  «  1. 


Substituting  (5)  into  (6)  gives  (dropping  the  Z  subscript 


on  n_) 


*  *  *  •  '  t  -  •  ■  t  • ■■  4V* 1  i  V  j..  i 


(7) 


% 


•’i* 


>1 


3 

|V 


,y» 


2/3  „  (Ze)2(4rrn)1/3 


which,  again,  is  the  plasma  parameter  relevant  to  an  OCP  with 


ions  of  charge  Ze 


The  canonical  expression  for  y  is  ('’"«n  by 


4ffnA: 


where  X_  is  the  Debye  distance 


13,14 


With  the  latter  two  expres¬ 


sions  we  find 


4wn(Ze) 


which  is  seen  to  be  the  Debye  distance  for  an  OCP  comprised  of 


ions  of  charge  Ze . 


These  relations  may  be  cast  in  terms  of  a  plasma  parameter 
more  common  to  studies  of  OCP.  It  is  given  by 


r  -  (Ze)  / akjjT 


(10) 


where  aJ  is  a  measure  of  the  mean  occupation  volume  per  ion. 


That  is, 


4  3  , 

j  ira  n  =  1 


(ID 


The  parallel  structure  of  y  and  F  is  evidenced  by  r^riting  (8) 


vm! 


Y 


(12a) 


(Ze)2/.\D 

ETT 


(12b) 

so  chat 


We  may  conclude  that  the  previously  stated  criterion  that 
separates  weakly  from  strongly  coupled  plasmas  may  also  be  given 
in  terms  of  r . 

Thus  we  find  that  in  the  limit  (4),  the  coupling  and 

parameters  of  a  two-component  plasma  of  electrons  and  ions  of 

charge  Ze  reduce  to  those  relevant  to  an  OCP  of  ions  of  charge 

Ze  in  a  negative  background.  We  may  conclude  that  thermodynamic 

properties  such  as  internal  energy  and  equation  of  state  for  a 

1  17 

two-component  laboratory  plasma  with  Z  »  Z  are  the  same  as 
those  of  an  OCP  comprised  of  the  same  species  of  ions. 

Numerical  work  of  Slattery,  Doolan  and  DeWitt  [9]  estab¬ 
lished  the  following  expression  for  internal  energy  U  for 

the  fluid  phase  of  an  OCP. 

fljpp  -  5r  +  br1/4  +  cr"1/4  +  d  +  Sr/ n  (13) 

B 

where  N  is  total  ion  number  and 

I  -  -0.398,  b  -  0.950,  c  =  0.190, 


d  -  -0.815.  5  =  0.010 


19 


The  Helmholtz  free  energy,  F,  may  then  be  obtained  through 
integration.  Namely, 
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fc.d  - 1 
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nn  =,  rr  flLCrJLl  +  3~ldr.  +  F<ri_> 

NkBT  J?1  NkBT  +  i  +  Uiqr 


(14) 


The  normalized  free  energy  F(T^)  /NkgT ,  ..itn  -  1  was  calculated 
employing  various  contributions  over  the  unit  interval.  There 
results  (in  the  limit  of  large  N)  [9] 

(15) 

Ijpf  -  ar  +  4(br1M  -  cr'1/4)  +  (d  +  3)  In  r  -  (a  +  4b  -  4c  +1.152) 

B 

With  this  value  of  free  energy  at  hand  an  equation  of  state  is 
obtained  from  the  thermodynamic  relation 


(16) 


To  perform  this  differentiation  we  first  rewrite  r  (10,12b)  in 
explicit  form 


(Ze)2  f4rTN]  1/3 


kBT  ^3  V 


(17) 


There  results 


(18) 


Differentiating  (15)  we  obtain 


1  3F  =  r- 

W^T  3  V  13 


»r'3/4  +  cr'5/4  +  (d-t-3 ) r ■ L ] 


(19) 


v'$| 


With  (16)  and  (13)  we  then  obtain 


PV 


l  +  |  ^  (5r  +  br1/4  +  cr‘I/4) 


(20) 


With  the  previously  stated  argument  we  may  conclude  that  (13) 

and  (20)  are  valid  energy  and  equation-of -state  formulas  for  a 

laboratory  plasma  comprised  of  ions  and  electrons  in  equilibrium 

at  a  given  temperature  and  obeying  the  constraint  (4) .  Note  in 

117  2 

particular  that  for  Z  =  5  (completely  ionized  boron)  Z  ‘  /Z 
0.26.  For  Z  » 13  (completely  ionized  aluminum)  this  ratio 
becomes  0.12.  Thus  the  proposed  equivalence  should  be  valid  for 
Z  >  5.  It  is  also  important  to  note  that  this  equivalence  is 
not  appropriate  to  processes  where  the  dynamics  of  electrons 
come  into  play,  such  as  conductivity  [12]. 


21 


ft# 


References 


1.  D.  M.  Heffernan  and  R.  L.  Liboff,  J.  Plasma  Phys.  27 ,  473 
(1982) . 

2.  K.  Brueckner  and  S.  Jorna,  Revs.  Mod.  Phys.  46 ,  325  (1975). 

3.  E.  E.  Salpeter  and  H.  M.  Van  Horn,  Astrophys.  J.  155 ,  183 
(1969) . 

4.  A  collection  of  papers  on  electromagnetic  propulsion  is 
given  in  IEEE  Trans,  on  Magnetics,  Mag-20 ,  March  (1984) . 

5.  S.  Ichimaru,  Revs .  Mod .  Phys .  54,  1017  (1982). 

6.  M.  Baus  and  J.  P.  Hansen,  Phys.  Rep.  59 ,  1  (1980). 

7.  G.  Kalman  and  P.  Carini,  Strongly  Coupled  Plasmas  (Plenum, 
New  York,  1978) . 

8.  E.  P.  Wigner,  Trans.  Faraday  Soc .  (London)  (1938). 

9.  W.  L.  Slattery,  G.  D.  Doolen  and  H.  E.  Dewitt,  Phys.  Rev. 
A21,  2087  (1980);  Phys.  Rev.  A26 ,  2255  (1982). 

10.  R.  L.  Liboff,  J.  Appl .  Phys.  56,  2530  (1984) 

11.  P.  M.  Platzman  and  P.  A.  Wolf,  Waves  and  Interactions  in 
Solid  State  Plasmas  (Academic,  New  York,  1973). 

12.  L.  Spitzer,  Jr.,  Physics  of  Fully  Ionized  Gases,  (Inter¬ 


science,  New  York,  1956). 


UNIFIED  EQUATION  OF  STATE  FOR  WEAKLY  AND  STRONGLY  COUPLED  PLASMAS 


Summary 


An  argument  is  presented  which  permits  a  laboratory  plasma 
of  arbitrary  ionization  to  be  viewed  as  a  one-component  plasma. 
In  this  equivalence,  effective  ion  charge  number  Z  is  dependent 
on  temperature  and  ion  number-density.  The  plasma  parameter  F 
thus  gains  additional  dependence  on  these  parameters  due  to  its 
dependence  on  Z.  Employing  previous  results  for  the  internal 
energy  of  such  fluids  gives  the  following  unified  equation  of 
state  for  a  laboratory  plasma  valid  over  the  interval, 


0  <  F  <  300. 


PV  -  Nk0T  + 


U*c(D0 


1  ? 


+  usc(r)0r  -  7 


u 


Here  we  have  written  U*^  and  ,  respectively,  for  weakly  and 
strongly-coupled  excess  energy  and  Q(x)  represents  the  unit 
step  function.  The  volume,  pressure,  temperature,  and  number 
of  ions  in  the  plasma  are  written  V,  P,  T  and  N,  respectively. 
A  numerical  plot  of  this  equation  as  a  function  of  temperature 
reflects  the  shell  structure  of  atoms  in  the  plasma. 
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INTRODUCTION 


The  study  of  strongly-coupled  plasmas  is  relevant  to  x-ray 

1  2 
lasing,  inertial-confinement  fusion  devices,  the  interior  of 

3 

certain  super-dense  stars,  and  in  plasma-driven  rail-gun 
devices.^  Studies  of  strongly-coupled  plasmas  have,  for  the 
most  part,  ’  ’  utilized  a  model  due  to  Wigner0  in  which  the 
plasma  is  viewed  as  ions  moving  in  a  uniform  charge-neutralizing 
background.  This  model  is  called  a  one-component  plasma  and 
carries  the  abbreviation  OCP. 

9 

In  an  earlier  work  a  means  of  constructing  the  equation 
of  state  for  a  laboratory  plasma  was  described  which  permitted 
use  of  OC  results. ^  However  the  consistency  of  this  analysis 
required  a  two-component  fluid  and  therefore  addressed  fully 
ionized  plasmas  only. 

In  the  present  work,  stemming  from  an  analysis  of 
Zel'dovich  and  Raizer,^  the  results  of  the  aforementioned 
study  are  extended  to  a  laboratory  plasma  of  arbitrary  ioniza¬ 
tion.  The  work  of  Zel ' dovich-Raizer  permits  the  plasma  to  be 
viewed  as  a  two-component  species  comprised  of  electrons  and 
ions  with  effective  charge-number,  Z.  The  parameter  Z  emerges 
as  an  implicit  function  of  temperature  T  and  ion  density,  n. 

Thus  the  equation  of  state  is  found  to  contain  additional 
temperature  and  density  dependence  in  Z  as  it  appears  in  the 
effective  plasma  parameter,  T. 
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// 


The  present  work  employs  f -dependent  expressions  for  the 
internal  energy  previously  obtained  for  strongly  coupled^ 

( 1  <  r  <  300)  and  weakly  coupled^  (F  <<  1)  OCP.  The  smooth 
connection  of  these  two  curves  near  T  ~  0.5  motivates  the 
extrapolation  of  the  strongly-coupled  form  to  the  domain 
0.5  <  T  <  1  This  connecting  segment  results  in  a  unified 
equation  of  state  valid  over  the  entire  domain,  0  <  T  <  300. 


Two  interesting  effects  run  through  the  analysis.  The 
first  of  these  pertains  to  the  form  of  the  effective  plasma 
parameter 

1/3 


-  n 

r  « 


(T.n) 


For  the  most  part,  previous  plasma  studies  did  not  view  Z  as 
temperature  dependent.  Under  such  circumstances,  the  plasma 
parameter  grows  large  with  decrease  in  temperature.  However, 
in  the  present  study  Z  is  seen  to  increase  monotonical ly  with 
T  (with  steps  at  ionization  shells)  and  we  find  an  overall 
decrease  of  T  with  T.  The  second  interesting  observation  is 
that  the  shell  structure  of  atoms  in  the  plasma  is  reflected 
in  the  equation  of  state. 

As  the  analysis  is  dependent  on  the  details  of  ionization 


energies,  a  specific  material  must  be  chosen  prior  to  numeri- 

13 


cal  work.  The  present  study  addresses  aluminum.  Thus  plots 
of  PV/NkgT  vs  T  as  well  as  "f  (where  P  is  pressure)  are  pre- 


15  20  -3 

sented  for  n  in  the  range  10  -  10  cm  v .  Plots  are  also 


included  of  Z  vs  T. 


25 


ANALYSIS 


Review 

Previous  Finding 

In  ref.  9  it  was  noted  that  for  a  two-component  plasma 
one  may  write  the  mean  two-particle  potential  energy  as 


<V>  s 


k 


4tt 

T 


1/3 


ni/Je‘ 

e 


-  (nfin) 


LQ. 


(1) 


where 
sity . 


ng  and  n  are,  respect 
With  conservation  of 

<V>  = 


ively,  el 
charge , 


ectron 
ng  =  2n, 


and  ion  number  den- 
(1)  becomes 


(2) 


providing 

Z2  »  Z1* 17  (3) 

We  recognize  (2)  to  be  the  average  interaction  potential  of  a 
OCP  comprised  of  ions  of  charge  Ze.  A  key  parameter  in  the  study 
of  strongly-coupled  plasmas  is  the  plasma  parameter,  which  with 
(2)  is  written 


r 


<v> 

V 


1/3 


(Ze) 


kBT 


(4) 


Thus  a  plasma  is  strongly  coupled  for  f  >  1  and  weakly  coupled 
for  r  <<  1 . ^ 

With  these  observations  and  employing  a  previously  obtained 
expression  for  the  internal  energy  of  an  OCP,^  an  equation  of 
state  was  obtained  relevant  to  the  domain  1  <  7  <  300. 
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Thermodynamic  Relations 


Prior  co  generalizing  these  results  to  a  laboratory  plasma 
of  arbitrary  ionization  we  present  a  brief  review  of  basic  ther¬ 
modynamic  relations  relevant  to  a  plasma. 

An  equation  of  state  for  a  plasma  is  related  to  internal 
energy  in  tl.e  following  manner  (see  Appendix  A). 

T* 


PV  ,  ir  .  i 
NFjT  3NkgT  l 


(5) 


where  we  have  set 


U*  -  U  -  |  NkgT 


(6) 


1 


Here  U  denotes  total  plasma  energy  and  U*  "excess"  internal 
energy.  Note  that  if  U*  is  known,  then  with  the  thermodynamic 
relation 


jHF 

av 


II  il 
ar  av 


i 


and  (4) ,  we  may  write 


?  PV 

r  a 

r 

F 

1  ^ 

1  Tr 

W  J 

where  F  is  the  Helmholtz  free  energy.  Integrating  the  preceding 
equation  at  constant  T  gives,  with  (5), 


fr  dr  (  u* 

J  r  ^ 


*  ' 
+  3 


F(r)  -  F(ro) 

N^x 


(7) 


Thus  knowledge  of  U*  as  a  function  of  r  for  a  plasma  gives  both 
the  equation  of  state  (5)  and  the  free  energy  (7). 


Effective  ionization 

As  noted  in  ref.  11,  a  plasma  in  equilibrium  at  a 

given  temperature  is  characterized  by  ions  of  effective  charge 
eZ,  where  e  is  (positive)  electron  charge.  Thus,  the  results 
of  ref.  9  carry  over  with  Z  replaced  by  an  effective  ion  charge 
number,  Z.  Specifically  one  may  write  (1)  with  Z  replaced  by 
Z  and  again  we  find  reduction  to  (2)  relevant  to  an  OCP,  in  the 
domain  ( 3 ) . 

The  ion-charge  number,  Z(T,n)  is  obtained  in  the  following 

manner.  First  the  discrete  ionization  energies,  I(Z),  as  a 

function  of  degree  of  ionization,  Z,  are  connected  by  line  seg- 

13 

ments,  such  as  depicted  in  Fig.  1  for  the  case  of  aluminum. 

This  gives  the  continuous  function  T(Z) .  The  function  Z(T,n)  is 
then  obtained  by  solving  the  following  implicit  equation,  which 
stems  from  the  Saha  equation  as  well  as  conservation  of  charge, 


Z 


at3/2 


exp 


I  Z  + 


i _ 

T 


\ 

1 

7 


(8) 
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where  A  =  6  x  10  cm  eV  and  T  is  in  eV.  In  the  present 
work  this  procedure  for  obtaining  Z(T,n)  was  carried  out  again 
for  the  case  of  aluminum.  Results  are  shown  in  Fig.  2  for  ion 
densities  (10l3-1020)  cm"3. 

The  plasma  parameter  7,  as  given  by  (4),  now  includes  addi¬ 
tional  temperature  and  density  dependence  through  Z(T,n).  That 
i  s  , 
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(9) 


r 


1/3 

[Z(T,n)e] 


2 


A  plot  of  this  parameter  vs  T  is  shown  in  Fig.  3. 


Equation  of  State 

Ordinarily  for  an  OOP,  U*  =  U*(F).  To  incorporate  our 
findings  as  described  above  in  an  equation  of  state  we  write 
U*  *  U*(7) ,  where  T  is  given  by  (9)  with  Z  determined  from  (8). 
Thus  (5)  becomes 


PV  ,  U*  ( n  , 

Nlc^T  3n1cbT 


(10) 
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In  the  weakly  coupled  domain  (F  <  0.5),  U  was  obtained  as  a 

8/2 

function  of  e  -  /3  F  .  Converting  these  results  to  a 
T-dependent  function  gives 


uSc(r> 

TnrBT 


a1r3/2  *  bj^r3  In  r  +  n^r3  +  djJ^lnr  +  (11a) 
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In  the  strongly-coupled  domain,  over  the  interval  1  <  F  <  300, 
we  have 


ir 

usc 


a2F 


+  b 


a/4 

2‘ 


+  c2r 


1/4 


+  d/ 


(lib) 


Constants  have  the  following  values 
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a1  -  -0.866  ,  bl  »  -1.125  ,  cl  -  -1.102 

dl  -  -2.923  ,  eL  -  0.243 

a2  -  -0.898  ,  b2  -  0.950  ,  c2  -  0.190  , 

d2  -  -0.815 

Inserting  these  results  into  (10)  gives  two  equations  of 
state  relevant,  respectively,  to  the  two  said  domains.  Numeri¬ 
cal  plots  (see  Fig.  4)  of  these  results  strongly  suggest  an 
interpolation  insert  over  the  interval  0.5  <  T  <  1,  given 
namely  by  the  strong  coupling  form  (11). 

With  this  interpretation  at  hand  one  is  able  to  write  a 
unified  equation  of  state  valid  over  the  whole  interval, 

0  <  f  <  300.  Namely, 


PV  - 


u*c<r>  a 


-  u*c(r>  e 


[r 


(12) 


where  7,  U*c  and  UgC  are  given  respectively,  by  (9)  and  (11), 
and  0(x)  is  written  for  the  unit  step  function  9(x)  ■  1,  x  >  0 
and  zero  elsewhere. 

In  application  of  (12)  one  should  recall  that  r  is 
implicitly  dependent  on  T  and  n  through  its  dependence  on  Z 
given  by  (8) .  A  numerical  plot  of  PV/NkgT  vs  T  obtained  from 
(12)  is  shown  in  Fig.  5.  Note  that  the  shell  structure  of 
atoms  in  the  plasma  is  reflected  in  this  equation  of  state. 
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CONCLUSIONS 


A  unified  equation  of  state  was  constructed  for  a  labora¬ 
tory  plasma  of  arbitrary  ionization.  This  was  achieved  by 
introducing  an  effective  ion  charge-number,  Z(T,n),  obtained 
from  the  equations  of  Saha  equilibrium.  This  form  of  Z  was 
then  used  to  obtain  a  generalized  plasma  parameter,  f(T,n). 
Employing  this  value  of  F  in  previously  obtained  expressions 
for  internal  energy,  and  effecting  an  interpolation  over  the 

internal  <  T  <  1,  gave  a  unified  equation  of  state  over  the 

—  —2  —1  17 

total  interval  0  <  r  <  300  for  effective  ionization  Z  »  Z 

It  is  important  to  note  that  the  equivalence  between  a  labora¬ 
tory  and  one-component  plasma  described  in  Section  A.l  is  not 
relevant  to  processes  where  the  dynamics  of  electrons  come  into 
play,  such  as  conductivity.^ 


-V'T-  - 

W! 
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Appendix  A 


In  this  appendix  we  wish  to  derive  the  relation  (5)  rele¬ 
vant  to  an  OCP.  First  we  recall  the  statistical  mechanical 

,  16 
relations 


W?  *  7  +  OTJf  C  u<r>S<r>4’r2dr 


( A  •  1 ) 


^3l-^rCu'(r)8(r)^3dr 


(A. 2) 


In  these  expressions,  g(r)  is  the  radial  distribution  function, 
and  u(r)  is  the  two-particle  interaction  potential.  Let  us 
assume  the  form 


u(r)  *  a  r 


(A. 3) 


where  a  and  s  are  constants.  Inserting  the  derivative 


u ’ (r)  -  - 


into  (A. 2)  gives 


PV 


1  +  C  u(r)g(r)^c:dr 


(A. 4) 


which  may  be  combined  with  (A.l)  to  give 


PV  s  U  _  3  I 


(A. 5) 


For  the  Coulomb  interaction,  s=l,  and  (A. 5)  returns  (5). 
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REVIEW  OF  ELECTRICAL  CONDUCTIVITY  TN  PLASMAS 


Summary 


A  review  of  studies  of  electrical  conductivity  in  a  plasma 
is  presented.  A  brief  description  of  domains  and  parameters  of 
plasma  physics  encountered  in  studies  of  conductivity  is 
included.  An  historical  account  of  conductivity  calculations 
is  given  and  it  is  concluded  that  one  such  study  is  most  rele¬ 
vant  to  conductivity  in  a  strongly-coupled  classical  plasma. 
Expressions  and  numerical  results  stemming  from  this  study  are 
included . 


Introduction 
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The  main  thrust  of  this  report  addresses  state-of-the-art 
studies  of  conductivity  in  a  classical  strongly-coupled  plasma. 

Preparatory  to  this  discussion,  in  part  II  of  the  report, 
plasma  domains  are  introduced  and  various  important  parameters 
are  discussed  relevant  to  the  natural  separation  of  these 
domains.  A  brief  description  of  the  fluid  picture  for  a  plasma 
is  included  as  well  as  a  discussion  of  domains  relevant  to  the 
case  where  a  magnetic  field  is  present. 

In  part  III  various  plasma  models  are  introduced  and  a 
discussion  of  Saha  equilibrium  is  included. 

Part  IV  presents  a  brief  historical  review  of  electrical 
conductivity  calculations  for  a  plasma. 

In  part  V  formulas  and  results  of  numerical  evaluation 
for  electrical  conductivity  are  presented.  Values  of  electri¬ 
cal  conductivity  as  a  function  of  the  plasma  and  compression 
parameters  are  given  in  a  table  at  the  conclusion  of  the  report. 

Plasma  Domains 
Parameters 

Physical  domains  of  a  plasma  divide  into  four  basic  areas^ 


as  deoicted  below. 


Weakly 

Coupled 

r  «  i 


Strongly 

Coupled 

r  >  1 


Degenerate 

rQ  ~  1 


The  parameters  in  the  above  inequalities  are  defined  as 
follows.  The  quantum  degeneracy  parameter  A  is  given  by 


Aq  -  n  Ad  u 

where  n  is  particle  number  density  and 

A  2  =■  „  h2  ( ' 

'  d  2:rml<gT 

is  the  thermal  deBroglie  wavelength. 

-1/3 

Thus,  in  the  quantum  domain  >  n  ,  and  the  deBroglie 
wavelength  is  of  the  order  of  interparticle  separation. 


W!\% 


iVt- 


Wi 


The  plasma  parameter 


rze)2 


* bT 


(3) 


(4ra3/3  -  n"1) 


is  a  measure  of  the  ratio  of  mean  two-particle  ion  potential  to 
particle  kinetic  energy.  Ion  density  is  written  n^ •  When 
f  «  1,  mean  kinetic  energy  dominates  potential  interaction  and 
the  plasma  is  weakly  coupled.  In  the  limit  F  >  1  mean  potential 
energy  begins  to  dominate  and  the  plasma  is  strongly  coupled. 


For  extreme  potential  dominance,  r  »  1,  and  one  expects  the 

2,3 


medium  to  undergo  a  phase  change. 


The  quantum  parameter  Fq  is  defined  as  follows' 


(4) 


TF 


where 


TF  "  £  _  2 

o^ne 


(5) 


is  the  Thomas  Fermi  length,^’3  E-  is  the  Fermi  energy  and  s  is 


the  dielectric  constant  of  the  medium.  In  the  weakly  coupled 


degenerate  domain  ■Q  «  1  and  densities  grow  so  large  that 


kinetic  energy  due  to  the  uncertainty  principle  dominates  poten- 

A 


tial  energy. 
Fluid  Picture 


The  plasma  frequency  is  given  by  (in  cgs) 


(6) 


Let  v  denote  the  dominant  collision  frequency  of  plasma  con¬ 
stituents.  When 


V  »  U, 

p 


(7) 


collisions  tend  to  diminish  collective  plasma  behavior  and  fluid 

7  8 

behavior  ensues.  ’  This  situation  is  depicted  below. 


Magnetic  Domains 

If  a  magnetic  field  B  permeates  the  plasma  then  additional 
criteria  come  into  play.  We  recall  that  a  free  particle  of 
charge  e  and  mass  m  undergoes  circular  motion  in  a  magnetic 
field  with  frequency  (cgs) 


n 


eB 

me 


(8) 


If  plasma  particles  suffer  collisions  at  frequency  v,  then  at 
v  »  SI,  the  circular  magnetic  motion  is  lost  to  collisions  and 
if  further  v  »  cu^  the  plasma  is  fluid  like. 

If  an  electric  field  E  is  also  present  then  with  Ohm's 


law  (cgs) 


2  >  V 


V  >>  Jj^  ,  -4  ! 

P  1 


Magnetic  j  |  Nonmagnetic 

JB/c  >  cE  i  |  JB/c  «  cE 


As  previously  noted,  a  plasma  with  a  magnetic  field  is 


fluid-like  for  v  »  w  ,,1.  If,  further,  JB/c  >  oE,  then  we  may 

D 


term  the  plasma  magnetofluid  dynamic.  Plasmas  with  v  >  ft  are 

Q 

.  ^  1  .  _  1 _ 1 _ _ _ *  _  /  WTTTN  \  J 


typically  termed  magnetohydrodynamic  (MHD) . 


An  assortment  of  waves  may  propagate  in  a  plasma  in  these 

....  , _ 10,11,12 


various  domains. 


Fluid  Picture  and  TransDort  Coefficients 


In  concluding  this  section  we  estimate  che  conditions 


under  which  a  plasma  is  fluid  like  with  respect  to  specific 
transport  processes.  In  making  this  estimate  we  recall  that 


H 


different  collision  processes  pertain  to  different  transport 
coefficients.  Thus,  for  example,  in  calculation  of  viscosity 

g 

in  a  plasma,  ion-ion  collisions  play  the  major  role.  For  this 
case  we  take  the  magnitude  of  the  Coulomb  ion-ion  cross  section 
to  be 

ou  =  t  [  (Ze)  2/kgT]  2  =  t  f2  a2  (10) 

This  gives  the  mean  free  path 

l  =  1 /nza i £  (11) 

and  collision  time 


r  -  1/ * kBT/M 


(12) 


where  M  and  n„  are  mass  and  number  density  of  the  ion  species. 
Combining  these  relations  and  recalling  (3,6)  we  obtain  (with 
n  -  Zn^) 


2  2  16  1  1 

)  1  T  mZ  75 


(13) 


We  may  conclude  that  for  calculation  of  viscosity,  a  fluid  pic¬ 
ture  is  appropriate  for 

r  >  21.4  (A/Z)L/3  (14) 


where  A  is  atomic  mass. 

For  evaluation  of  electrical  conductivity,  electron-ion 
8 

collisions  dominate  and  we  write 


and 


<kBT)2 
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(15) 


%  * 


when  acorns  are  stripped  of  all  electrons  leaving  bare  nuclei 
and  electrons  in  the  plasma. 


Saha  Equilibrium 

The  equilibrium  state  of  a  laboratory  plasma  is  described 


by  the  Saha  equation 


15,16,17 


2uz  exp(-AEz/kBT) 


*Z-1 


u 


Z-l 


(19) 


n  \7 
e  d 


* 
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» 
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V 
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We  have  written  nz  for  the  density  of  ions  of  charge  eZ  and 
AEz  is  the  ground  state  ionization  energy  for  the  transition 
Z-l  -*  Z.  The  electronic  partition  function  uz  is  given  by 


1 

UZ  =  1 

4  i-0 


(i) 


'  -E 


(i) 


°z  exP!"¥?r 

'  D 


(20) 


where  gz*^  and 


are,  respectively,  statistical  weight  and 
.  th 


excitation  energy  of  the  i  excited  state  of  the  Z  ion.  The 
series  (20)  terminates  at  the  value  i*  corresponding  to  orbits 
whose  radii  are  comparable  with  mean  distance  between  constitu¬ 
ents  of  the  plasma.  It  has  been  demonstrated  by  Zel'dovich  and 
Raizer^^  that  for  a  majority  of  ions,  the  ground-state  term  in 
(20)  dominates  over  the  remaining  terms.  Thus  (20)  reduces  to 


‘Z-l 


2g<°> 

«Z-l 


exp ( -aE7/k3T) 

- n — 

n  \  . 
e  d 


(21) 


'•8 


W  it 


I 


(16) 


In 

'  target 


1 


this  latter  expression  we  have  recalled  that  it  is  the 
ion  density  which  is  relevant.  There  results 


2  2  16  M  Z3 

up  T  “  T  5  3 


(17) 


Thus,  £or  a  fluid  picture  to  be  relevant  for  electrical  conduc¬ 
tivity  one  must  have 


r  >  21.4  ZA1/3  (18) 

The  inequalities  (13)  and  (17)  follow  from  the  fluid-picture 
criterion  (7) . 


Plasma  Models 
PC?  and  Other  Models 

A  model  often  employed  to  describe  a  strongly  coupled 
plasma  is  the  so-called  one-component  plasma,  which  carries 

the  abbreviation,  OCP .  This  model  is  comprised  of  ions  moving 

1  ^ 

in  a  charge-neutral  background.  It  was  conceived  by  Wigner 
to  study  phase  change  to  the  solid  state.  An  extensive  review 
of  the  physics  of  an  OCP  has  been  given  by  Ichimaru.^  This 
model  is  often  employed  in  the  study  of  strongly  coupled 
plasmas . 

The  simplest  laboratory  plasma  model  is  that  of  a  fully- 
ion ized  or  two-component  plasma  comprised  of  electrons  and  ions 
of  charge  Ze .  This  situation  is  approached  in  the  limit  of 
very  high  temperature  (kgT  »  I^,  I^  3  atomic  ionization  energy) 

49 
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This  equation  indicates  that  in  an  actual  laboratory  plasma 
at  finite  temperature,  all  ionic  species  of  atoms  are  present 
at  varying  relative  densities. 

Historical  Review 

Earliest  expressions  for  electrical  conductivity  in  a 

90  91 

plasma  are  due  to  Spitzer  and  Harm“  and  Braginskii.  These 
results  are  relevant  to  a  fully  ionized  nondegenerate  weakly- 
coupled  plasma. 

These  findings  were  improved  upon  by  Rogers,  DeWitt  and 
Boercker"  in  a  study  of  electrical  conductivity  in  a  partially 

ionized  non-degenerate  plasma.  In  this  calculation  the  authors 

93  1 A 

employed  a  Chapman-Enskog-  expansion  (with  Sonine  poly- 

25 

nomials)  of  the  Wi 1 1 iams -DeWi tt  kinetic  equation.  Scatter¬ 
ing  from  electronic  shells  of  ions  and  neutrals  was  treated 

7  fi 

using  effective  interaction  potentials  developed  by  Rogers. 

We  note  that  the  Williams-DeWitt  equation  is  the  standard 
3oltzmann  equation  with  quantum  corrections  entering  through 
the  scattering  cross  section.  It  is  sometimes  called  a 
3alescu-Lenard  equation  due  to  shielding  structure  included  in 
the  interaction  potential. 

n  *7 

In  the  work  of  Boercker  et  al.Ll  a  plasma  collision  fre- 

28  ">9 

quencv  is  developed  which  correctly  reduces  to  the  Ziman*" 

2  5 

and  3alescu-Lenard^  results  in  appropriate  limits.  The  Ziman 
formulation  considers  a  strongly-coupled  degenerate  plasma  at 


so 


0  3  X  . 
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The  work  of  ref.  27  is  a  generalization  of  Boercker’s 

3  1 

in  which  a  c-matrix  formulation  of  the  Kubo  formula  is 

employed  to  calculate  electrical  conductivity. 

Working  with  the  Uehl ing-Uhlenbeck  quantum  generalized 

33 

Boltzmann  equation  Lampe ,  in  an  early  study,  employs  a 
Chapman-Enskog  expansion  about  the  Fermi-Dirac  distribution 
to  calculate  thermal  conductivity,  <,  which  may  be  related  to 
mobility,  u ,  through  the  Einstein  relation 


Here  we  have  written  cv  for  specific  heat  per  particle  in  the 
plasma . 

Lee  and  More0  employ  a  Krook-Bhatnager-Gross  equation 
to  calculate  transport  coefficients  in  a  plasma.  This  technique 
depends  on  knowledge  of  the  relaxation  time  which  the  authors 
calculate  working  with  the  Fermi-Dirac  equilibrium  distribution. 


Electrical  Conductivity 
Strongly  Coupled  Classical  Plasma 

The  most  relevant  expression  for  electrical  conductivity 
in  a  strongly  coupled  quas iclass ical  plasma  among  those 

described  in  the  preceding  section,  appears  to  be  that  obtained 

27 

by  Boercker  et  al. 


It  is  given  by 


where  r.  the  plasma  parameter,  is  given  by  (3).  With 
3  =  (koT)’1-,  the  collision  frequency  v  is  given  by 


(  13/2  ,  dk  k3  e"3f'2k‘'/8m  u  .  (k)v  (k) 

mn  ,31  c 30  e  l  e  i 

"  *  tl  irrml  i,  - rrFToi - 


(23) 


^See(k)Sii(k)  '  Sei(k)Sie(k)1 


A  closely  allied  expression  for  v  was  obtained  by  Baus  et  al . 

In  (23),  Sa^(k)  is  the  (dimensionless)  dynamic  structure  factor 

relevant  to  the  two  components,  a  and  b.  The  electron-ion 

Fourier  transformed  potential  is  v  .  (k)  and  the  Fourier  trans- 

e  i 

formed  generalized  potential  is  uQ^(k).  These  potentials  have 
dimensions  of  energy  x  volume.  The  transform  ugi(k)  satisfies 


the  relation 


uab<k)  *  -cab(k)/3 


where  c  K(k)  is  the  Fourier  transform  of  the  direct  correlation 

3  D 

T  C  17 

function.  ’  The  factor  \(D  in  (22)  is  a  correction  factor 
to  account  for  the  fact  that  v  as  given  by  (23)  is  equivalent 
to  a  s ingle-Sonine  polynomial  approximation. 

A  rough  estimate  of  \(f)  is  obtained  in  the  weak-coupl ing 


limit  which  gives 


.  93  oj  (3"/2) 1/2 


"A 


(2A1 


,-3/2, 


(Recall  that  in  cgs  units  a  has  the  dimensions  of  frequency.) 
The  function  A  is  given  by 


A  = 


-L  -  dke-af’2k2/8%  Su(k) 


- - 7  J 

4  it  e  0 


£(k,o) | 


(25) 


In  the  Debye-Huckel  limit  one  obtains 


A  =  e1/ct  E]_  (1/a)  -  j  el/2a  Ej.  (l/2a) 


■j  (In  a  -  y  -  In  2)  =  -y  (In  a  -1.27) 


(26) 


where 


8m 

a  =  - , 

SfrkD 

E^(x)  is  the  exponential  integral  and  y  is  Euler's  constant. 
Numerical  Results 

Numerical  evaluation  of  A  given  by  (25)  was  performed 

using  direct  correlation  functions  and  structure  factors  from 

solution  of  the  hypernetted  chain  equation  (HNC).  With  these 

values  of  A  inserted  into  (22,23),  numerical  integration  for 

27 

a/w0  was  performed  at  two  values  of  the  compression 
parameter^ ' ® 


r _  s  a /a 
s  o 


where  aQ  is  the  Bohr  r^uxus.  Values  of  o/jj  stemming  from  the 
Debye-Huckel  result  (26)  were  also  obtained.  These  values  are 


SA 


_ 


listed  in  Table  I. 
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Conclusions 


ft 
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A  review  of  studies  of  electrical  conductivity  in  a  plasma 
has  been  presented.  The  report  began  with  a  description  of 
domains  and  parameters  of  plasma  physics  common  to  such  studies. 
An  historical  account  of  calculations  led  to  the  conclusion 
that  the  formulas  most  germain  to  conductivity  in  a  classical 
strongly  coupled  plasma  were  those  obtained  by  Boercker  et  al 
These  findings  were  reported  in  section  V  concluding  with  a 
list  of  numerical  results  given  in  Table  1. 
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NONLINEAR  ELECTRICAL  CONDUCTIVITY  FOR  STRONGLY  COUPLED  PLASMAS 


Summary 


A  nonlinear  analysis  of  electrical  conductivity  in  a  plasma 
is  given  stemming  from  the  Uehling-Uhlenbeck  equation.  Aniso¬ 
tropy  due  to  an  applied  electric  field  is  incorporated  through 
a  Legendre  polynomial  expansion  of  the  distribution  function. 

The  plasma  is  comprised  of  ions,  electrons  and  a  neutral  compo¬ 
nent.  The  electron-ion  interaction  is  described  by  a  shielded 
Debye  potential  at  high  energy  and  a  cut-off  Coulomb  potential 
at  low  energy.  A  nonlinear  equation  for  the  distribution  func¬ 
tion  is  solved  and  yields 


£sl(x) 


1  +  Be 


A  (x) 


for  the  symmetric  part  of  the  solution.  Nondimens ional  energy 
is  x,  B  is  a  normalization  constant  and  A(x)  is  an  explicit 
integral  dependent  on  the  electric  field  and  specifics  of  the 
interaction.  Resulting  nondimensional  conductivity  a,  is  given 
by 


1  r2  3/2  VZ+1) 
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dx 


where  Z  is  effective  ionization,  is  the  ratio  of  charge  to 
total  heavy-particle  density,  Q  is  dimensionless,  weighted 
cross  section  and  Aq  and  F^  are  quantum  and  plasma  parameters 
respectively.  Application  is  made  to  an  aluminum  plasma  and 
plots  of  conductivity  vs  electric  field  are  obtained.  These 

plots  exhibit  three  distinct  regions.  With  increase  in  field 
strength  these  are:  Ohmic,  Coulomb-dominated  and  neutral- 
dominated  . 
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Introduction 

Previous  kinetic  studies  of  electrical  conductivity  in 
plasmas  have,  for  the  most  part,  involved  linear  theory.  Thus, 
for  example,  approaches  based  on  the  Kubo  ’  formula  involve 
linear  response  theory;  those  based  on  a  Chapman-Enskog  expan¬ 
sion  employ  a  linearized  collision  integral  in  the  Boltzmann 

3  4  S 

equation;  ’  whereas  others  start  with  a  linear  equation  such 

as  the  linear  Krook  equation,^  or  the  linearized  Fokker-Planck 
equation.^  As  a  consequence,  such  analyses  can  obtain  expres¬ 
sions  for  current  which,  at  best,  are  linearly  dependent  on 
electric  field. 

The  present  non-linear  analysis  begins  with  the  Uehling- 

Q 

Uhlenbeck  quasi-clcssical  kinetic  equation.  A  Lorentz  expan- 
a 

sion  of  the  distribution  function  in  terms  of  Legendre  Poly¬ 
nomials  ^ ^ ^  is  employed  to  account  for  anisotropy  in  the 
distribution  function  due  to  the  presence  of  an  electric  field. 
The  analysis  addresses  a  three-component  plasma  comprised  of 
electrons,  heavy  ions  of  arbitrary  ionization  Z,  and  neutrals. 
These  latter  two  components  are  assumed  to  be  in  equilibrium. 
The  interaction  between  electrons  and  ions  is  described  by  a 
shielded  Coulomb  potential.  It  is  found  that  the  presence  of 
the  neutral  component  inhibits  over  population  of  the  tail  of 
the  distribution  function,  better  known  as  the  'runaway 
effect '  .  ^ 


Two  key  parameters  which  enter  the  analysis  are  the  Debye 
plasma  parameter  and  quantum  degeneracy  parameter,  Aq.^ 

Thus  when  >  1  the  plasma  is  termed  strongly  coupled  and  when 
Ap  >  1  it  is  degenerate.  The  present  ..uric  bridges  the  quasi- 
classical  (Aq  =  1)  and  classical  domains  (Aq  «  1).  Further¬ 
more,  apart  from  the  assumption  of  a  specific  interaction  poten¬ 
tial,  the  analysis  should  prove  valid  in  the  strongly-coupled 
domain  as  well. 

A  nonlinear  equation  is  obtained  for  the  symmetric  part 
of  the  distribution  and  solved  exactly.  In  the  limit  of  zero 
electric  field  the  total  distribution  reduces  to  the  Fermi  Dirac 
distribution.  For  smal 1 -electric  field  it  gives  the  displaced 
Fermi  Dirac  distribution.  In  the  classical  domain  it  returns 
the  Druyvesteyn  distribution . ^  ^ 

Integrating  the  total  distribution  function  gives  a  closed 
expression  for  conductivity  as  a  function  of  electric  field. 
Numerical  integration  at  constant  heavy-particle  temperature, 
density  and  ionization,  yields  conductivity  vs  electric  field 
for  various  values  of  ion  density.  For  reasonable  charge  den¬ 
sities  three  distinct  regions  are  evident.  At  low  field  values 
Ohm's  law  is  obeyed.  At  intermediate  fields  the  conductivity 
rises  due  to  effects  of  the  Coulomb  cross  section.  At  high 
fields  the  hard  core  potential  of  the  neutral  component  domi¬ 
nates  yielding  the  familiar  E  fall  off  of  conductivity.^ 


Starting  equation 

Our  starting  equation  is  the  Uehling-Uhlenbeck  quasi- 

8 

classical  kinetic  equation,  appropriate  to  a  plasma  comprised 
of  electrons,  ions  and  neutrals.  It  is  assumed  that  ions  and 
neutrals  are  of  the  same  mass,  M,  and  are  both  in  equilibrium 
at  the  temperature  T.  These  equilibrium  distributions  are 


given  by 


VV  "  acW 


W  "  a»W 


where 


w 


.  I  M 


3/2  -Mvl‘/2kgT 


and  a~n  and  a»Tn  are  ion  and  neutral  number  densities  respec- 
u  o  No 

tively.  With  aN  + ac  *  1,  nQ  represents  total  heavy-particle 
number  density. 

With  these  assumptions  our  starting  kinetic  equation  for 
the  electron  distribution  function,  f(x,v,t)  is  given  by 

4l+v.l=+??*4s-  J^CE . F_>  +  J„C£.F„)  ( 


where  E  is  applied  electr-'~  field  and  m  is  electron  mass.  The 


.  17. 


collision  integrals  are  given  by  (in  Boltzmann  notation  ) 


J.(f,F)  -  [[f'(l-Cf)F[  -  f Cl-Cf ,)F1]a.gdndv1 


(4a) 


where 

C  =  (2irfi/m)  ^ ,  g  -  f  v  -  Y]_  I  (4b) 

and  F  denotes  either  F^  or  F^  and  a ^  is  the  corresponding  elec¬ 
tron  heavy-particle  cross  section.  Inserting  (1)  into  the  RHS 
of  (3)  gives 

3(f)  2  +  JN  -  J[f,(l-Cf)F(^1  -  fd-cf'JFQj^JogdQdVj^  (5) 

Here  o  represents  the  composite  cross  section 

a  *  aCaC  +  aN°N  (6) 


Combining  (5)  and  (3)  gives  our  starting  kinetic  equation  for 
the  subsequent  analysis.  We  will  apply  this  equation  to  an 
equilibrium  homogeneous  plasma  for  purposes  of  calculating  elec¬ 
trical  conductivity.  There  results 


eE 

"of 


3f 

3v 


3(f) 


with  J ( f )  given  by  (5). 


(7) 


Lorentz  Expansion 

To  account  for  anisotropy  in  the  distribution  f(v)  we 
expand  it  in  Legendre  polynomials  as  follows. 

oo 

f(v)  =  Z  P.(y)f.(v)  (8) 

i=o 

where 

u  =  v  •  E  (9) 
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and  hatted  variables  are  unit  vectors.  Substituting  (8)  into 
the  LHS  of  (7)  and  keeping  £*1  terms  gives 


df,  2f ,  df 

-  3v  T  “3F  ~  +pi“3T 


Substituting  the  series  (8)  into  the  RHS  of  (7)  we  obtain 


(1C) 


J(fQ  +  yf1)  -  J(fQ)  +  AJ(fo,f1) 


(11a) 


where 


AJ(f0,f1)  -  r { C  y • f ^ (1-Cf^-^uf 1 )  -  f;(5uf,)]F^ 


O  lov^n;jl01 


(lib) 


-  [uf1(l-Cf;-£y’fi)  -  fo(CM,f[)]F01}agdndv1 

[Note  that  =*  f^(v)  whereas  FQ1  *  FQ(v1).]  With 

2  .  o 


taken  as  a  parameter  of  smallness  we  find  (see  appendix  A) 

f(H  -  F01  *  0<E>  • 


v'  ■  v  +  0(  -: ) 


Keeping  terms  of  leading  order,  (lib)  becomes 


AJ(f0,fi)  -  n  /  (u'-u)cgdft 


With 


Q(v)  =  {(1-cos  3)  j(v,e)dn 


(12) 


(13) 


the  preceding  equation  becomes,  to  this  same  order  (see 
Appendix  3) 


_ 


AJ(f0,f1) 


-"of^vQ 


(14) 


Combining  (7,  10,  11  and  14)  and  equating  coefficients  of  like 
Legendre-polynomial  degree  gives 


eE 

35 


dfl  .  2fl 


J(f0) 


eE  o 


¥  "3v  -  -"oVS 

We  note  that  (15a)  may  be  rewritten 


(15a) 


(15b) 


35  ^  Hv  (v2fl)  "  J(fo} 


(16) 


Integrating  v  in  this  equation  over  a  sphere  of  radius  v  we  find 


v2fL(v)  -  I ( v) 


(17) 


where 


l(v)  =  j{/V  [f;(l-Cf0)F(J1  -  fo(l-?f;)F01]ogvZdv}dfidflvdv]L  (18) 


and 


dv  *  v4dvdQ. 


Evaluation  of  I(v) 

As  is  shown  in  Appendix  C,  I(v)  may  be  reduced  to  the  form 
(to  leading  order  in  s) 


I  (v) 


f\  r+iv 

J  Hv 


fo(l-Sf^)F01cgv“dv| dfidfivdvL 

j 


(19) 


where 


Av  -  v,  •  (g  -  g ' ) 


(20) 


We  concentrate  on  the  integrand  of  (19) 


X  I  fo(l-5£;)F0llJgv' 


(21) 


With  the  observations 


we  expand 


v  <  v  <  v  +  Av,  av  »  0(e) 


-  d£o(v> 

£0  ■  £0Cv)  *  £0<v)  +  <v_'')  '37~" 


df  (v) 

f;  -  f0(v,)  a  -Hhr-  + 


and  obtain,  to  leading  order  in  e, 


X  3  YLf0(v)  l  l-CfQ(v)  ]  +  ”■  2 


df  (v)  „  df  (v) 

*  Wv>  Hv 


(22) 


where  i  (v ,  v ,  8  ,  v^) are  known  functions  (see  Appendix  D.) 
Writing  (18)  to  leading  order  in  e  we  find 


1  *  Ilfo(1“Cfo)  +  l2  "3v"  +  1 3  ^ f o  Uv" 


(23) 


where 


r  r v+iv 

:j  5  JJv  *i 


dvd^v.dr»dv^ 


Following  a  technique  introduced  by  Davydov^  we  concentrate  on 

V 


r  r  V+AV  „  „  "5 

J  Jv  (v-v)Fq1(v1)c(v, 9)v3dvdflvdQdv1 


(25) 


Integrating  over  v  and  noting  the  equali  ty 


*  ^  A  , 

j  V,  v,  dv 


J  -1  -1  a~l  I  J  a-l 


f  dv. 


^where  T  is  the  unit  tensor)  we  find 


1 2  *  -j  /  Vj2  (1-cos  9)F01(v1)a(v,  0)v3dJ2vdndy1 


With  the  distribution  (2)  we  note  that 


f F0l(vl)  7  Mvl  d-l  ■  7  nokBT 


which  gives 


-  n  k-T 

4itv3  -2-L.  q(v) 


(26) 


It  is  evident  from  (18)  that  I(v)  depends  on  E  only  through 

its  dependence  on  fQ.  However  we  see  that  1^  as  given  by  (23,24) 

do  not  depend  on  f  .  With  this  observation  we  evaluate  I.  from 

o  j 

17)  by  setting  E  =  0 .  We  choose  fQ(v,E)  at  E  ■  0  to  be  the  Fermi 
Dirac  distribution 


f0(v.o) 


5  +  Bq  exp (mv^ /2ksT) 


(27) 


^hich  has  the  prooerties 


f)8 


Furthermore,  f  (v,E)  is  normalized  by 


\  f  4^v^dv  =  n 
o  0 


(29) 


where  n  is  electron  density.  With  this  distribution  at  hand  we 


obtain 


_  mv  T  T  r\ 

1  FTT  X2  *  X3  “  0 


Inserting  these  values  into  (23)  gives  the  result 

4:rn  v3Q(v)  [\  ,  „  df 

I(v)  =■  - g -  !nvf0(v)  (l-CfQ(v))  +  kgT  (v) 

With  I(v)  thus  evaluated  we  return  to  (17). 


(30) 


(31) 


Equation  for  f(v) 

Combining  (15b)  and  (17)  gives 


4*  feEi2  v  dfo 

T  tTj  "HT'  i(v) 


(32) 


Inserting  (31)  into  this  equation  and  introducing  the  dimension¬ 
less  energy, 


.  mv 

X  =  9U  T 


(33) 


gives,  after  some  reduction,  the  nonlinear  equation 


SR 


*1*1*1  'WWKVWM 


Integrating  (34)  we  find 


First  note  that  (36)  and  (37)  give 


df . 


[x+s(x)J  +  xf0(l-€f0) 


(34) 


c  .  I  M 

s  =  Z  5 


eE 


•s  ? 


(35) 


f0oo 


1  +  BeA<^ 


(36) 


A(x)  =  f 


x  x  '  dx  1 


+  s(x') 


(37) 


dfo 

“35T 


x+s  (‘1'4fo)fo 


(38) 


which  is  seen  to  reduce  to  (28)  in  the  limit  s *  0 . 

With  (38)  and  recalling  (8)  and  (15b) ,  we  obtain  the 
total  distribution  function 

f ( x )  -  fQ(x) ( l-u[6es(x) ] 1/2  fl-5f0(x)]}  (39) 

It  is  convenient  at  this  point  to  introduce  the  dimensionless 
distribution^ 


fSL(x)  =  Sf0(*> 


1  +  3e 


A  ( x ) 


(40) 


where  3  e  B/$.  This  distribution  has  the  normalization 
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M  "  J 


f°°  Tqr  (x) 1  ^dx 
'o 


(41) 


where 


A 


Q 


n  A 


3 

d 


and  is  the  thermal  deBroglie  wavelength 


(42a) 


.  2  =  2*fi~ 

Ad  =  mkgT 


(42b) 


Note  that  in  the  classical  limit, ^  C  -*■  0  [see  (4b)  ]  and  f^L(x) 

goes  to  the  Druyvesceyn  distribution.^  For  low  electric  field, 

from  (35,37)  we  see  that  s=0,  A(x)  =x  and  (38,39)  return  the 

19 

displaced  Fermi-Dirac  distribution. 


Conductivity 

Calculation  of  Q  integrals 

Employing  the  distribution  (39)  gives  the  current  density 
(after  an  integration  by  parts) 


J 


r e]  2  V 

H  no 


d 

cfx 


,  * 

'Q- 


TSL(x)dx 


(43) 


Evaluation  of  the  integral  in  this  expression  demands  knowledge 
of  Q(x)  as  given  by  (13).  Inserting  the  composite  cross  section 
(6)  into  (13)  we  find 


Q  *  j  (l  -  cos  3)  [  +  a^j^jdft 


or,  equivalently, 


Q  =  Qc  +  Qn 


(44a) 


(44b) 


71 


where  and  correspond  Co  Che  charge  and  neutral  components 

2( 

of  (44a),  respectively.  For  we  employ  the  results  of  Liboff 
for  the  high-energy  domain  and  those  of  Chapman  and  Cowling^ 
relevant  to  the  low-energy  domain.  In  the  Tirst  of  these  cal¬ 
culations  the  electron-ion  interaction  potential  was  taken  to  be 


’(r)  *  exp(-r/XD) 


(45a) 


where  the  Debye  distance  A-  is  given  by 


4ffa~n  e' 
L  o 


Z(Z+1) 


(45b) 


In  the  low-energy  calculation  the  potential  was  taken  to  be 


V(r) 


r  <  A, 


(46) 


and  zero  for  r  >  A^.  Related  cross  sections  are  as  follows. 
We  introduce  the  dimensionless  cross  section 


Q  i  7T  A2  Q 


(47) 


where 


.  .  Ze 
"  - 


is  the  so-called  "distance  of  closest  approach."  Liboff's 


result  may  then  be  written 


-§  h  % *  * 


J  =  In  : 


(48) 


m 


•  V*. 


where  y  is  Euler's  constant  and 


r  =  — 

D  '  XD 

is  the  conventional  plasma  parameter.  Chapman-Cowl ing ' s  result 
appears  as 


Qc,c-c  *  jp  ln 


1  + 


2x 


(49) 


A  smooth  interpolation  formula  bridging  the  two  resu  Its  (48,49) 
is  given  by 


•  ,  aC 

k  “7 


See  Fig.  1.  For  x  >  rQ,  (50)  gives  Q 
x  <  rD.  it  gives  Qc  =  Qc>c.c. 


1  + 

2 

+  3 

1  -  exp 

r  4  ^ 

X 

'TJ 

n 

_ 

_ 

(50) 


^  whereas  for 


Lastly  for  QN  we  set 
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(aa0) 


JN 


(51) 


where  a  is  a  pure  number  between  one  and  six  and  a  is  the  Bohr 

o 

radius.  For  example,  for  Al ,  a  =  2.4.  With  (51)  we  find 

'  2 


QN  "  aN 


aa. 


^  j 


(52) 


Substituting  (50)  and  v52)  into  (&4b)  gives  the  desired  com- 


With  (47),  che  expression  (35)  may  be  rewritten 


\  2 


where  the  characteristic  electric  field 


!bk  V4  V 

Eo  ■  Vir  - 5 - 


(53) 


(54) 


Furthermore  we  note  that  the  parameters  and  may  be  con- 

2  2 

veniently  written  in  terms  of  the  Bohr  radius,  a„  =  fi  /me  and 
J  o 

2 

Rydberg,  R  =  e  /2a  . 


rD  ^  ^ac23(Z+D(n0a03)f^]3 

D 


(„)3  2,21  .  3  I  :  3 

’  =  [  o  °  J  C<bTJ 


(55a) 


(55b) 


These  forms  give 


n  a„ 
o  o 


3  Z(Z+1) 


2rr  a 


(56a) 


C  1 D 


7  Z4(Z*1)2 


with  these  expressions  at  hand  Eq  as  given  by  (54)  may  be 
written 


(56b) 


2 

where  e/a J- 
o 


31/2  f m'| 1  / 2  Z5/3(Z+H1/3  A 


e/a  ~ 
o 


Q 

7271 


(57) 


5.14  x  10’  V/cm. 


Evaluation  of  Conductivity 
Rewriting  (43)  in  the  form 


J  *  o  E 


and  setting 


o  *  c/u> 


(58) 


(59) 


where 


0  4Tra-n  Ze' 
l  Co 


we  obtain 


,  fn)3/2  a-(Z+l)i/z  ,  f  ) 

12  CN  x  I  j 

T  T  - \ — r -  J  t  ct  (x)  ~XZ  “>  dx 

J  AQrD  o  SL  dx  [QJ 


(60) 


where  fgL(x)  is  given  by  (40)  and  Q  by  (46b)  .  (50)  and  (52) . 

Note  that  a  as  given  by  (60)  is  a  function  of  E,  Z,  a^ ,  nQ  and  T 

where,  we  recall,  T  is  heavy-particle  temperature.  We  apply 

20 

this  formula  to  an  aluminum  plasma  with  »  10  and  k„T  »  5eV 

O  D 

for  which  Z  =  2.5.  With  these  values,  (57)  gives 
Eq  *  8.97  x  10^  V/cm.  Taking  a  =»  2.4  in  (51)  permits  evaluation 
of  o  as  given  by  (60)  with  a^  a  free  variable.  Numerical  plots 
of  a  vs  (E/Eq)  for  these  values  at  a^  =  0.1,  0.5  and  0.9  are 
shown  in  Fig.  2. 


V 


Three  regions  of  behavior  are  apparent  from  these  curves. 

For  low  electric  field  the  curves  are  seen  to  obey  Ohm's  law. 

At  E  =  0 . 1Eq  nonlinear  effects  due  to  the  Coulomb  cross  section 

come  into  play  which  are  evidently  more  pronounced  for  higher 

ion  density,  a^n0.  At  E  =  0.3Eo  the  Coulomb  cross  section  gives 

-1/2 

way  to  the  cross  section  of  the  neutrals  resulting  in  a  E 
fall  off  of  conductivity. 


Conclusions 

Starting  with  the  Uehl ing-Uhlenbeck  quasi-classical  kinetic 
equation  we  have  addressed  the  problem  of  nonlinear  electrical 
conductivity  in  a  plasma  comprised  of  electrons,  ions  and 
neutrals.  A  Legendre  polynomial  expansion  was  introduced  to 
account  for  anisotropy  of  the  electron  distribution  function 
due  to  the  applied  electric  field.  An  interpolation  was  intro¬ 
duced  to  bridge  the  low  and  high  energy  components  of  the 
electron-ion  cross  section.  An  exact  solution  of  the  resulting 
equation  for  the  electron  distribution  was  obtained  which  was 
found  to  go  to  consistent  limits  in  the  classical  and  quantum 
limits.  With  this  distribution  at  hand  an  explicit  formula  for 


electrical  conductivity  was  obtained  in  terms  of  E,  Z,  ar ,  n 
and  T.  Numerical  plots  of  conductivity  vs  electric  field  for 
3  specific  case  of  an  aluminum  plasma  at  varying  vaiues  of  a^ 
were  obtained.  These  curves  revealed  three  distinct  domains  of 
behavior  which  were  associated  with  various  components  of  the 
cross  section. 


7b 


Appendix  A 


In  this  appendix  we  wish  to  justify  the  relations  pre¬ 
ceding  (12) .  With 


2  _  m 
e  "  M+m 

the  center  of  :cass  velocity  is  given  by 

G  =*  v  +  (l-e^)v^ 

With  relative  velocities  defined  by 

§  -  y  -  ?i 
g*  -  Y'  -  vj^ 

we  find 

V2  -  V2  +  2(l-t2)G  •  (§'-  g) 

v’2  -  v,2  -  2s2G  •  (g1-  g) 

2 

Dropping  terms  of  order  z  then  gives 

Yi  *  Cg '  -  g) 

v '  =  v  +  — - ^ - - — 


(Al) 


(A2) 


(A3. a) 
(  A3  .  b) 


( A4 . a) 
(A4.b) 


( A5 .  a) 


vj^  =*  v1  (A3.b) 

•7 

The  latter  of  these  relations  gives,  to  0(e~),  the  desired 
result 

F  '  »  F 

*01  ‘01 
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For  electron  temperature  Te  equal  to  or  in  excess  of  ion 
temperature  T,  we  write 


I  kBT  -  <7  Mv2>  S  <7  mv2>  *  7  kBT« 


or  equivalently,  <x>  >  1.  In  addition  this  relation  gives 


(A6) 


(A7) 


This  result  together  with  (A3. a)  permits  (A5.a)  to  be  written 


v’  =  v  + yL  •  (£'- g) 


(A8) 


which  to  0(e)  returns  the  relation 


v  =  v 


(A9) 
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TO** 


>  r*V 


Appendix  B 


Appendix  C 


In  this  appendix,  we  wish  to  derive  (19).  Recalling  (18) 
we  write 

=  /  |  r  U^l-U^F^  -  fo(l-i£;)F01]ogv2dvjd!Idnvdv1  (Cl) 

Consider  the  transformation  which  exchanges  primed  with  unprimed 
variables.  We  find 


=  lY-yJ 

-  lY'-YjJ  s  g- 

(C2a) 

dg  dy '  dy £ 

■  dg  '  dy  dv^ 

(C2b) 

<?(g) 

*  o(g’) 

(C2c) 

where  (C2c)  holds  due  to  the  invariance  of  the  inverse  cross- 
section.^  With  (C2)  at  hand,  we  obtain 

(C3) 

/  ^  fo(1^fo)F0l^YdQdv1  =  [[  ^  fo(l-Cfo)F01agdvdfldv1 

|  v  |  <V  j  V  '  I  Cv 

From  (A8)  |v'|  <  v  and  |v|  <  v  +  v^  •  (g-g')  are  equivalent  to 
order  e.  Combining  this  fact  with  (Cl)  and  (C3)  returns  the 
desired  result,  (19). 


Appendix  D 


In  this  appendix  the  functions  4^  of  (22)  are  found.  With 
the  observations 


v  <  v  <  v  +  Av  (DI) 

and  that  Av  as  given  by  (20)  is  of  order  e,  we  find 


g  “  V  +  0(s) 

(D2a) 

v4  -  vL  +  0(e) 

(D2b) 

v-v  -  0(e) 

(D2c) 

v'  -  v  =»  0(e) 

(D2d) 

a(g,g' )  -  a(v,g' )  +  0(e) 

(D2e) 

Inserting  these  expressions  into  (21)  and  expanding  fQ  and  f^ 
about  v,  we  obtain 

X  i  f0(l-£fo)FQ1agv2  -  fo(v)(l-Cfo(v))F01(v1)a(vI£' )v3 

df  „  o 

+  (v-v)  (v)(l-Cfo(v))F01(v1)a(v>g')vJ 


(D3) 


df. 


Cf0(v)(v'-v)  (v)F01(v1)a(v,g’ )v 


-3 


+  A(g’  ,vIv,y1)fo(v)  (l-^fQCv))Fo(v1)  +  0(e“) 


where  the  function  A(g  ,v,y,v^)  is  of  order  e.  Comparison  of 
the  above  expression  with  (22)  serves  to  identify  the  ^ 
functions . 
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